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1 Introduction 

Accurate modeling and simulation of turbulent flow is a topic of intense ongoing 
research. Two main approaches are identified, differing by the amount of detail 
that is included in the physical and numerical description. Direct numerical 
simulation (DNS) aims to calculate the full, unsteady solution to the governing 
Navier-Stokes equations. While accurate in principle, DNS is severely restricted 
by limitations in spatial resolution. This situation summons alternative simu- 
lation approaches that are aimed at capturing the primary features of the flow 
above a certain length-scale only. A prominent example of this is the large-eddy 
simulation (LES) strategy in which a smoothing of the flow features and a cor- 
responding reduction in the flow complexity is introduced by spatial filtering, at 
the expense of introducing a 'subgrid' closure problem. 

We consider so-called Leray regularization of the convective contributions 
Q] . This gives rise to a subgrid parameterization which involves both explicit 
filtering and (approximate) inversion. The Leray model also arises from the a- 
modeling strategy derived via Kelvin's circulation theorem [| . We study the 
dynamics associated with the Leray model in a turbulent mixing layer and com- 
pare predictions with filtered DNS results and findings due to dynamic (mixed) 
models || . In particular, the kinetic energy, momentum thickness and energy- 
spectra are analyzed, establishing favorable performance of the Leray model and 
robustness at arbitrarily high Reynolds number. This is unique for a similarity- 
type model that does not contain an explicit eddy-viscosity term. 

We provide the basic Leray formulation in section |^ together with the nu- 
merical inversion of the filter. Application to turbulent mixing is presented in 
section and concluding remarks are collected in section EL 
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2 Basic Leray formulation 

In the filtering approach the evolution of the filtered solution {ilj, p} is governed 
by the spatially smoothed Navier-Stokes equations. Filtering the nonlinear terms 
gives rise to the turbulent stress tensor Here m — L(u.i) is 

the filtered velocity field in the Xi direction, with L denoting the linear filtering 
operation which we assume to have a (formal or approximate) inversion L . 
Likewise, p = L(p) is the filtered pressure. Expressing Ty in terms of the filtered 
velocity is the basic closure problem in LES. 

Leray regularization provides an intuitively appealing method for modeling 
Tij. In this formulation the convective fluxes are replaced by UjdjUi, i.e., the 
solution u is convected with a smoothed velocity u. Consequently, the nonlinear 
effects are reduced by an amount governed by the smoothing properties of L. 
For commuting filters L the governing equations in the Leray formulation can 
be written as 

djUj = djUj = ; d t Ui + UjdjUi + dip — -j^-djjUi = (1) 

tie 

Uniqueness and regularity of the solution to these equations have been estab- 
lished rigorously pj . The Leray formulation contains the unfiltered Navier- 
Stokes equations as a limiting case. The unfiltered solution can readily be elim- 
inated, giving rise to a closed formulation for {Hi,p}: 

d t Ui + dj (HjUi) + dip - + dj m ij = ( 2 ) 

where the asymmetric, filtered similarity- type Leray model mfj arises: 

771^ = L^UjL^ 1 ^)^ - u-jlli (3) 

In LES, one commonly adopts compact support filters with filter-kernel G. 
In one spatial dimension such filtering can be expressed as 

/•A/2 

u{x) = L(u) = / G(z)u(x + z) dz (4) 

J-A/2 



where A is the filter-width. In actual simulations the resolved fields are known 

lm=0- 



only at a set of grid points {a; m }^ =0 . Numerical filtering corresponds to kernels 



G(z) = '£a j d(z-z j ) ; \ Zj \ < A/2 (5) 

We consider three-point filters with ciq = 1 — a, a\ = a_i = a/2 and zq = 0, 
z\ = —z-i = A/2. In addition, we use a = 1/3 which corresponds to Simpson 
quadrature of the top- hat filter. The application of L~ l to a general discrete 
solution {u(x m )} can be specified using discrete Fourier transformation as H] 



L-\u m ) = £ ( a - 1 + VT ^ f\l-2a)-^u m+rj/2 (6) 



j=-n 
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where the subgrid resolution r = A/h (with h the grid-spacing) is assumed to be 
even. An accurate and efficient inversion can be obtained with only a few terms, 
recovering the original signal to within machine accuracy. At fixed A, variation 
of the subgrid resolution r allows an independent control over flow-smoothing 
and numerical representation H . 



3 Shear layers at arbitrary Reynolds number 

The turbulent mixing layer [|| is simulated in a volume L 3 at various Re. At a 
modest Re = 50, an assessment of the quality of Leray modeling is obtained by 
comparing with filtered DNS data and dynamic models. Moreover, we consider 
this flow at high Re, adopting a fourth order accurate spatial discretization. 




Figure 1: Momentum thickness 9 (a) and streamwise kinetic energy spectrum 
A at t — 75 (b): Leray (32 3 : solid, 64 3 : dotted), dynamic model (32 3 : dashed), 
dynamic mixed model (32 3 : dash-dotted), filtered DNS (o). LES at A = L/16. 

Visualization of the DNS, obtained on 192 3 cells, demonstrates a well devel- 
oped flow beyond t = 40. In Fig. |l|(a) the evolution of the momentum thickness 
shows Leray predictions to compare well with filtered DNS data and with dy- 
namic (mixed) models at 32 3 . This is improved using 64 3 points. In Fig. |l|(b) 
the streamwise kinetic energy spectrum at t = 75 on 32 3 shows improved cap- 
turing of the large and intermediate scales, compared to dynamic models. The 
dissipation of the smaller scales is significantly improved on 64 3 . The Leray 
model shows both forward and backward scattering of energy [^j . 

A particularly appealing property of Leray modeling is the robustness at very 
high and even infinite Reynolds number, cf. Fig. || (a,b). Although comparison 
with filtered DNS data is impossible, we observed that the smoothed flow dy- 
namics is properly captured and a nearly grid-independent solution is obtained 
as r — A/h > 4 || . The kinetic energy is likely to be slightly over-predicted. 
At high Re the spectrum displays a region with k~ 5 ^ 3 behavior. 
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Figure 2: Kinetic energy E (a) and streamwise kinetic energy spectrum A at 
t = 75 (b): Leray model at Re = 50 (solid), Re = 500 (dashed), Re = 5000 
(dash-dotted). LES at 32 3 and A = L/16. The dotted line represents k~ 5 / 3 . 

4 Concluding remarks 

The Leray model was found to predict the momentum thickness properly while 
exhibiting both forward and backward transfer of energy. Further analysis shows 
reliable levels of turbulence intensities and correct behavior of kinetic energy. 
The Leray model has a tendency to underestimate dissipation. The computa- 
tional overhead associated with the Leray model is lower than that of dynamic 
(mixed) models and no introduction of ad hoc parameters is required. The reg- 
ularized dynamics shows an appealing robustness at high Re. 
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